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Abstract: The High Energy cosmic Radiation Detection (HERD) facility is a space mission designed for detecting 


cosmic ray (CR) electrons, 7-rays up to tens of TeV and CR nuclei from proton to iron up to several PeV. The main 


instrument of HERD is a 3-D imaging calorimeter (CALO) composed of nearly ten thousand cubic LYSO crystals. A 
large dynamic range of single HERD CALO Cell (HCC) is necessary to achieve HERD’s PeV observation objectives, 
which means that the response of HCC should maintain a good linearity from minimum ionizing particle (MIP) 


calibration to PeV shower maximum. In order to study the linearity of HCC over such a large energy range, a beam 


test has been implemented at the E2 and E3 beam lines of BEPC. High intensity pulsed electron beam provided by 


E2 line are used for producing high energy density within HCC; 1*/proton provided by E3 line are used for HCC 


calibration. The results show that no saturation effect occurs and the linearity of HCC is better than 10% from 30 
MeV (1 MIP) to 1.1x10? TeV (energy density is 93 TeV/cm*), which can meet the requirement mentioned above. 
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1 Introduction 


HERD is one of the space astronomy payloads of the 
cosmic light house program onboard the China’s Space 
Station. The main scientific objectives of HERD are in- 
direct dark matter search, precise CR spectrum and com- 
position measurements up to the knee energy and high 
energy y-ray monitoring and survey [1]. To achieve these 
scientific objectives, a 3-D imaging and five-side active 
calorimeter is designed to perform high energy resolution 
and high statistics measurements of the CR nuclei, elec- 
trons and positrons, and 4-rays in space. HERD CALO 
consists of 21x21x21 cells corresponding to 55 radiation 
lengths (Xo) and 3 nuclear interaction lengths (Az) longi- 
tudinal, separately. Each HCC is made of a 3x3x3 cm? 
cubic LYSO crystal and two spiral WLSFs (Saint Gobain 
BCF914A) as the fluorescence readout. All the WLSF sig- 
nals are collected by two image intensified CCDs (ICCD) 
2. 


LYSO (Lutetium Yttrium Orthosilicate with Cerium 


HERD, Calorimeter, LYSO, Linearity, Beam test 


doping) crystal is one of the best candidates for CALO 
material. It has the advantages of high density, high 
light output, short decay time, and application of this 
crystal in space would be quite easier thanks to the non- 
hygroscopicity and small temperature coefficient. The 
general properties of LYSO are described in Ref. 
Theoretically, the fluorescence intensity of LYSO crystal 
should be linear with the energy deposition induced by 
incident particles. However, saturation effect may occur 
at high ionization density of the crystal. The linearity 
of LYSO crystal light output should be studied and un- 
derstood before application. Results from Monte Carlo 
simulation show that a dynamical range of 2x10° is re- 
quired for HCC [I]: A minimum detectable signal down 
to 30 MeV is required for calibration using MIP; mean- 
while, the maximum energy deposition induced by PeV 
proton shower is up to 60 TeV (energy density is about 
8 TeV/cm?, Fig. |1). 
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Fig. 1: Energy deposition in 1x1x1 cm? LYSO cubes 
with 1 PeV proton incidence. The whole calorimeter con- 
sists of 100x100x100 cubes corresponding to 5 A; deep, 
thus the shower maximum produced by 1 PeV proton 
can be contained within the calorimeter for those early 
showered events. The data is derived from Monte Carlo 
simulation using Fluka 2011 with DPMJETS3 [6]. 


Numerous studies of LYSO linearity have been per- 
formed using low energy ^-rays or electrons, the verified 
ranges are 511 keV—1333 keV [7], 65 MeV—145 MeV 
and 100 MeV—1.5 GeV [9]. There are few calorimeters 
of space missions designed to measure PeV particles di- 
rectly, thus little attention has been paid to the linearity 
of inorganic crystals over a large energy (or energy den- 
sity) range. 

A study of linearity for LSO (Lutetium Orthosilicate 
with Cerium doping) crystal using high intensity pulsed 
y-rays is described in Ref. [10], which demonstrates that 
the linear response of the crystal ranges from 5.1x 10/5 
MeV /(cm?.s)—1.9x10!? MeV/(cm?-s). It is reasonable 
to expect the same linearity of HCC which is made of 
similar material at that high energy density. However, 
the verified energy density of this study is too high and 
it does not cover the required range of HCC. 

We developed a method of using high intensity pulsed 
electron beam and 7* /proton beam to measure the lin- 
earity of HCC. In this paper, the method and the exper- 
imental setup are presented in Section [2] Data analysis 
and results are presented in Section[3|and the conclusion 
is briefly shown in Section [4] 


2 Experimental setup 


2.1 Test Beam 


It is well known that ground accelerators can only 
produce primary beam particles up to several TeV, 
thus it is feasible to achieve an energy density up to 8 
TeV/cm? by using a bunch of particles passing through 
HCC and calibrate it by using single beam particle, sep- 


arately. The measurement is carried out by using the 
electron beam provided by BEPC E2 Line [I1]. The 
beam properties are showed as follows: 


e Beam energy: 2.5 GeV 

e Beam intensity: from 10? to 10'° electrons/pulse 
e Beam frequency: 12.5 Hz 

e Beam size: 4 x 4 cm? 

e Beam pulse width: 20 ps 


The collimation size is 1 cm which is smaller than crys- 
tal size (3 cm) in order to reduce the geometrical non- 
uniformities. Energy deposition of one electron with 2.5 
GeV kinetic energy in HCC is about 240 MeV (Sec- 
tion [3.3]. HCC is thick enough (2.6X,) for 2.5 GeV 
electrons to produce showers, the part of HCC at the 
shower center has the highest energy density of about 
20 MeV/cm?. Thus an adequate beam intensity up to 
4x10? electrons/pulse (or energy deposition in HCC up 
to 96 TeV) is necessary to get a maximum energy density 
of 8 TeV/cm?. 

BEPC E3 line is a mixed beam of secondary particles 
includes y, e*, 7*, proton and w*. m*/proton mixed 
beam with the momentum of 600 MeV/c to 1.2 GeV/c 
are chosen by beam line magnets for HCC calibration. 


2.2 The HCC readout 


PMTs instead of ICCDs are used during the test to 
avoid possible non-linearity of the ICCD device. Two 
WLSFs coupled to a cubic LYSO crystal (from Suzhou 
Jtcrystal Co., Ltd.) are read out by two PMTs (XP2262) 
respectively. Neutral density filters (NDF) with differ- 
ent transmittances are attached to PMTs for detecting 
different energy ranges. The NDFs and the PMT are 
all sealed in a black plastic cylinder box with a circular 
aperture of 0.4 mm on top for the WLSF to plug in, as 
shown in Fig. 
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2.3 Experimental setup at E2 line 


Three detectors are used to calibrate the beam inten- 


sity (Fig. Bp: 


e Faraday Cup (FC): FC is the only device which 
can measure the absolute beam intensity directly. 
However, FC is not sensitive enough, it is designed 
to monitor the beam with intensity 10$ elec- 
trons/pulse; 


e Ionization Chamber (IC): designed to monitor the 
beam intensity from 10^ to 107 electrons/pulse; 


e Thin Plastic Scintillator (PS): placed in front of 
the HCC and attached with 2 PMTs for readout 
purpose. The effective detecting range is from 10? 
to 10° electrons/pulse. 


With cross calibration between the three detectors, a 
wide beam intensity from 10? to 10? electrons/pulse can 
then be monitored and precisely measured. 
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Fig. 3: Schematic diagram of experimental setup for 


linearity measurement at E2 line. 


A lead wall is built for shielding the detectors from the 
contamination caused by secondaries from the interac- 
tion of beam particles with surrounding materials. Sig- 
nals from FC and IC are recorded by a 12-bit peak sens- 
ing ADC (Mod. V785N), while signals from 4 PMTs are 
readout by a 4-channel digitizer (Mod. DT5751). Both 
DAQ modules share a common trigger signal from the 
accelerator. 


2.4 Experimental setup at E3 line 
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Fig. 4: Schematic diagram of experimental setup for 
HCC calibration at E3 line, the NDFs are removed and 
the high voltages for PMTs are raised to detect the MIP 
signal which is much weaker than electron signal. 


In the configuration at E3 line, a time of flight sys- 
tem (TOF) consists of two plastic scintillators are used 
for 7* /proton separation as shown in Fig.|4| The trigger 
is the coincidence of two PMTs attached to the plastic 
scintillators. DT5751 which can work at a sampling rate 
of 1 GS/s provides a high time resolution to distinguish 
proton signals and 7* signals from TOF (Section [8.3]. 


3 Data analysis 


3.1 Calibration of beam intensity monitoring 
devices 


A cross-calibration between FC, IC and PS is per- 
formed as shown in Fig. 


HCC 


lonization 
Plastic C i 


10? 10° 10* 10° 10° 10’ 10° 10? 


Electrons/Pulse 


1 10 


Fig. 5: Principle of beam intensity calibration. Overlaps 
between the detecting ranges of FC, IC and PS are large 
enough for a cross-calibration. 


3.1.1 FC Calibration 
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Fig. 6: Relationship between injected electron number 
and FC output, the x axis is ADC channel; the y axis is 
number of electrons. 


Fig. [6] shows good linearity in ADC's counting range, the 
calibration factor is obtained by fitting a straight line to 
data points: 


1 FC Channel = 2.58 x 10* +185 electrons (1) 


3.1.2 IC and PS Calibration 


The IC and PS detectors are calibrated by study- 
ing the relationship between FC output and IC output 
(Fig. (7), IC output and PS output (Fig. [). 
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Fig. 7: Profile histogram of IC calibration for displaying 
the relationship between FC output and IC output. 


ADC channels of IC from 200 to 750 are chosen for fit- 
ting, because when ADC(IC)>750, rare events would 
lead large fluctuations, and when ADC(IC) «200, signals 
would be too close to the baseline. 
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Fig. 8: PS calibration, the x axis represents PS out- 
put (Digitizer channel); the y axis is IC output (ADC 
channel), the histogram shows good linearity. 


'The calibration factors of IC and PS are obtained by 
fitting to the data: 


1 IC Channel = 3697 + 26 electrons (2) 


1 FC Channel = 21.25 +0.15 electrons (3) 


3.2 NDF Calibration 


A measurement of transmittance of NDFs was per- 
formed before the beam test, we used an optical fiber 
irradiated by a LED as the light source and a PMT for 
readout. The NDF was installed on a linear stage which 
could provide high precision linear movement. The orig- 
inal signal and attenuated signal could be obtained by 
moving the NDF to the position between the light source 
and the PMT. Fig. [9] shows the calibration result of one 
NDF. 
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Fig. 9: Calibration result of one NDF with 0.6 optical 
density(O.D., Transmittance = 107 ?-P-). 


Table 1: 


Calibration results of the NDFs used in the test. 


No. Attached PMT Diameter (mm) Optical density Nominal transmittance Calibrated transmittance 
1 L 504 0.6 0.251 0.2326 
2 L 504 2 0.01 0.0118 
3 H 504 il 0.1 0.0828 
4 H 504 3 0.001 0.0010 


The total transmittance of attached NDFs are 
2.8x 107° and 8.3x 107" (Table[1) for the two PMTs (de- 
fined as PMT.L, PMT.H), thus the intensity of light 
reached to the windows of PMT_L is about 34 times 
higher than that reached to PM'T H. The dynamic range 
of the HCC readout system is then extended naturally 
by combining the two PMTs. 


3.3 HCC Calibration using E3 beam 


Most of E3 beam particles are protons and 7*, TOF 
data which clearly illustrates the difference between them 
in Fig. [10] is needed to separate 7* from proton to get a 
clean spectrum as shown in Fig. 
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Fig. 10: Time of flight distribution measured by TOF 
system for E3 beam at 800 MeV/c. The unmarked peaks 
are caused by accidental coincidences. 
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Fig. 11: Use the TOF data to separate 7* from protons 
and use landau convoluted Gaussian function to fit at 
peak. The incident particle is identified as a 7 * when 10 
ns« TOF «15 ns or a proton when 16 nsc TOF «21 ns. 


Energy deposition in HCC for different particles is 
then calculated by Monte Carlo simulation using Geant 
4.9.6 [12]. The simulation environment is built based on 
the experimental setup. The results are shown in Fig. 
and Table [2] 

Energy deposition is nearly constant for 7+ in the 
case of normal incident when its momentum 400 
MeV/c; for 2.5 GeV electrons, energy deposition is much 
higher because electromagnetic showers are initiated in- 
side the crystal. 


S 914 EERH 2.5GeV electron 
E EN |2GeV/c x 

N x L] 1.2GeV/c proton 
E 01 E 0.8GeV/c proton 
3 [| 0.6GeV/c proton 


0.08 


95 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 0.26 0.28 0.3 


Energy Deposition in HCC(GeV) 
Fig. 12: Energy deposition in HCC for zt, protons and 
electrons derived from Mote Carlo simulation data. 


3.4 Results and Discussions 


'The results of linearity at different beam intensities 
as well as correlation between output of PMT_L and 
PMT_H are illustrated in Fig. 

To calculate the equivalent HCC output for proton and 
mt, the NDF?’ attenuate effect and the change of high 
voltage for PMTs should be taken into account: 


HCC,+ proton — Any Ar Di proton (4) 


Where Apy represents the high voltage correction 
factor which is defined as the ratio of two gains of PMTs 
used by E2 test and E3 test; Ap represents the total 
transmittance of NDFs used in E2 test; S, proton repre- 
sents the PMT output given by E3 test (TableB). Fig. [14] 
shows the profile histogram which combines all the re- 
sults for different intensity ranges as well as 7* and pro- 
tons. 
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Fig. 13: Linearity of HCC at different beam intensities (a,b,c) and correlation between output of PMT_L and PMT.H 
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Table 2: Summary of simulation results and experimental results of beam test at E3. 
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Fig. 14: Linearity of HCC (a) and deviations from the fit- 
ting line (b) for all ranges. The magnets at E2 cannot be 
precisely controlled to achieve an intensity lower than 200 
electrons/pulse, thus a huge gap without measurement 
points is located in between electrons and 7* /protons. 


The maximum deviation is defined as linearity of 
HCC. Fig.{14|shows the linearity is better than 10% when 
energy deposition (energy density) in HCC are from 205 
GeV (17 GeV /cm? ) to 1.1x 10? TeV (93 TeV /em?) which 
covers the requirement. The zt and proton signals agree 
well with the fitting function. 

There are some factors that may affect the measure- 
ment result of linearity: 


a) Non-linearity of PMTs: Non-linearity of two PMTs 
was studied before the beam test. The details of 


4079 428.2 
11644.6 783.1 
7796.4 630.9 
5715.4 564.8 


measurement are described in [I3]. Result shows 
the non-linearity of PMT is less than 296 when the 
output signal is not higher than 2 V. NDFs with 
high attenuation are used in the beam test to en- 
sure the output signals are lower than 2 V. 


Stability of PMTs: Before data acquisition, the 
high voltage of PMTs are usually switched on at 
least two hours in advance, the response of PMTs 
gradually become stable in this period. The tem- 
perature variation may also cause a drift. To min- 
imize such effect, one full data acquisition is com- 
pleted within 3 hours and the data acquisition is 
repeated several times. The temperature variation 
during the beam test is smaller than 2°C. There- 
fore, this residual influence is negligible comparing 
to the calibration error. 


Backgrounds at E2 line: The backgrounds are sec- 
ondary particles from the interaction of beam elec- 
trons with shielding material. Most of the back- 
grounds are X-rays and y-rays. One part of them 
are absorbed by the LYSO crystal. But the ratio 
of signal to background is lower than 2% (Fig. 
which do not have significant impacts on the mea- 
surement. 


Signal 


rrr 


HCC output from PMT_H (Digitizer channel) 
TTET [ 
co 
D 
oa 
: = 
o 
£ 
3 
+ Q 


& 


ZW a m L a u.s. acl 
105 10° 107 
Beam intensity (Electrons/pulse) 


Fig. 15: The background data is obtained by mov- 
ing HCC 3 cm away from the beam line. In this 
case, only secondary particles can hit on HCC. 
Both beam signal and background have a linear 
relationship with the beam intensity, the ratio of 
signal to background is about 1.7%. 


d) Calibration Error: Because the FC has a large 
size, perfect shielding is nearly impossible. Pri- 
mary electrons may hit on FC directly and lead 
to an overestimation of beam intensity. The large 
fluctuation of FC response can also have an influ- 
ence on the accuracy of IC and PS calibration. 


4 Conclusion 


In this paper, the linearity of LYSO crystal for HCC 
is studied by measuring the relationship between HCC 
output and electron beam intensity given by faraday 
cup, ionization chamber and plastic scintillator. The 
linearity is better than 1096 from 205 GeV to 1.1x10? 
TeV (or from 17 GeV/cm? to 93 TeV /cm?) which meets 
the requirement of PeV hadron detection. The linear 


range can be extrapolated to 7* (MIP, 30MeV) and pro- 
tons. No obvious saturation effect of HCC is observed in 
this range. The uncertainty of linearity measurement is 
mainly dominated by the measurement system and beam 
monitoring system. 

From this work, we get a preliminary conclusion that the 
LYSO has a good linearity and can be used in HERD 
CALO for detecting high energy particles. A beam test 
at CERN SPS, which can supply electrons and protons 
with energies up to 300 GeV, is planned to verify main 
performances of HERD CALO. 
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